Structure of pair winds from compact objects with application to 
emission from hot bare strange stars 
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ABSTRACT 

We consider a stationary, spherically outflowing wind consisting of electron-positron pairs 
and photons. We do not assume thermal equilibrium, and include the two-body processes that 
occur in such a wind: M0ller and Bhaba scattering of pairs, Compton scattering, two-photon pair 
annihilation, and two-photon pair production, together with their radiative three-body variants: 
bremsstrahlung, double Compton scattering, and three-photon pair annihilation, with their in- 
verse processes. In the concrete example described here, the wind injection source is a hot, bare, 
strange star. Such stars are thought to be powerful sources of pairs created by the Coulomb bar- 
rier at the quark surface. We present a new, finite-difference scheme for solving the relativistic 
kinetic Boltzmann equations for pairs and photons. Using this method we study the kinetics of 
the wind particles and the emerging emission for total luminosities oi L = 10'^'* — 10^^ ergs s~^ 
(the upper limit being set, at the moment, by computational limitations). We find the rates 
of particle number and energy outflows, outflow velocities, number densities, energy spectra, 
and other parameters for both photons and pairs as functions of the distance. We find that for 
L > 2 X 10'^^ ergs s^^, photons dominate the emerging emission. For all values of L the number 
rate of emerging pairs is bounded: ^ N^'^^ — 10"*^ s~^. As L increases from ^ 10^"' to 10*^ 
ergs s~^, the mean energy of emergent photons decreases from ~ 400 — 500 keV to 40 keV, as the 
spectrum changes in shape from that of a wide annihilation line to nearly a blackbody spectrum 
with a high energy (> 100 keV) tail. These results are pertinent to the deduction of the outside 
appearance of hot bare strange stars, which might help discern them from neutron stars. 
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There is now compelling evidence that electron- 
positron (e^) pairs form and flow away from many 
compact astronomical objects identified as neu- 
tron stars and black holes. Among these objects 
are pulsars (Sturrock 1971; Ruderman & Suther- 
land 1975, Arons 1981; Usov & Melrose 1996; Bar- 
ing & Harding 2001), accretion-disc coronas of the 
Galactic X-ray binaries (White & Lightman 1989; 
Sunyaev et al. 1992), active galactic nuclei (Jel- 
ley 1966; Herterich 1974; Begelman, Blandford, 
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& Rees 1984; Lightman & Zdziarski 1987; Sikora 
1994; Blandford & Levinson 1995; Reynolds et 
al. 1996; Yamasaki, Takahara, & Kusunose 1999), 
cosmological 7-ray bursters (Eichler et al. 1989; 
Paczyhski 1990; Usov 1992), etc. Several pair- 
production mechanisms are thought to operate 
in these objects. In pulsar magnetospheres, it 
is supposedly effected via conversion of single 7- 
photons into pairs in a strong magnetic field 
{'jB e+e^i?), while in galactic nuclei, the 
pairs are created mostly in photon-photon colli- 
sions (77 —^ e~^e^). Thermal neutrinos and an- 
tineutrinos radiated by a very hot compact ob- 
ject may be absorbed near its surface and pro- 
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very young (age < 10 s), compact objects with 
mean internal temperatures of ~ 10^^ K, the pair 
luminosity produced via this mechanism may be 
as high as ^ lO^'^ ergs s~^ (e.g., Eichlcr ct al. 
1989; Haensel, Paczynski, & Amsterdamski 1991), 
which is a typical luminosity of cosmological 7- 
ray bursters. Such pair winds are, therefore, dis- 
cussed as possible sources of cosmological 7-ray 
bursts (for a review, see Piran 2000). Recently 
it was shown that strange stars, which arc made 
entirely of deconfined quarks (e.g., Witten 1984; 
Alcock, Farhi, & Olinto 1986a; Haensel, Zdunik, 
& Schacffer 1986), may also be powerful sources 
of pairs (Usov 1998, 2001a). In this case, the 
pairs are created in an extremely strong electric 
field at the quark surface. The thermal luminos- 
ity of a bare strange star in pairs depends on 
the surface temperature Tg and may be as high as 
^ 10^° — lO^'^ ergs s~^ at the moment of forma- 
tion when Ts may be up to ~ 10^^ K (see Fig. 1). 
The luminosity of a young bare strange star in 
pairs may remain high enough (> 10^^ ergs s~"^) 
for ~ 10^ yr (Page & Usov 2002). 

The properties of pair plasmas have been stud- 
ied extensively, but most of the initial attention 
focused on plasmas in thermal equilibrium (e.g., 
Bisnovatyi-Kogan, Zeldovich, & Sunyaev 1971; 
Lightman 1981; Svensson 1982; Guilbert & Step- 
ney 1985; Zdziarski 1985 and references therein). 
It is conceivable that such equilibrated plasmas 
exist in some objects. Very powerful pair winds 
are among these objects. For an energy injection 
rate of E > 10^'^ ergs s^^, the pair density near a 
compact object of radius ~ 10^ cm is very high, 
and the outflowing pairs and photons are nearly in 
thermal equilibrium almost up to the wind photo- 
sphere (e.g., Paczynski 1990). The outflow process 
of such a wind may be described fairly well by 
relativistic hydrodynamics (Paczynski 1986, 1990; 
Goodman 1986; Grimsrud & Wasserman 1998; 
Iwamoto & Takahara 2002). The emerging emis- 
sion consists mostly of photons, so ~ E. The 
photon spectrum is roughly a blackbody with a 
temperature of ~ 10i°(E/105° ergs s"!)!/"* K. The 
emerging luminosity in pairs is very small, 
ie = E - L-y ~ 10~''L^ < L^. All this applies 
roughly down to E ~ 10"^^ _ iQiS 

ergs s . For 

E < 10^^ ergs s~^, which is the region we explore, 
the thermalization time for the pairs and pho- 
tons is longer than the escape time, and pairs and 



photons are not in thermal equilibrium (see be- 
low). Recently, a brief account of the emerging 
emission from such a pair wind has been given by 
Aksenov, Milgrom, & Usov (2003). In the present 
paper, we study numerically the kinetics of 
pairs and photons in a stationary pair wind with 
energy injection rates E < 10"*^ ergs s^^, and find 
the structure of the wind and its emerging emis- 
sion in pairs and photons. We assume that the 
outflowing wind consist of pairs and photons 
only. For the sake of concreteness, we take, in the 
present study, the input pair injection at the inner 
surface to correspond to a hot bare strange star. 
We use the Boltzmann equations for both pairs 
and photons (kinetic theory approach). The prin- 
cipal advantage of this approach compared to the 
Monte Carlo method is that it gives better photon 
statistics. Besides, the kinetic theory approach 
may still be used when the optical thickness is 
^ 1, where the Monte Carlo method is ineffective. 
Some examples of the kinetic theory approach for 
investigations of non-equilibrium pair plasmas can 
be found in Fabian et al. (1986), Ghisellini (1987), 
Svensson (1987), and Coppi (1992). The remain- 
der of the paper is organized as follows. In § 2 
we formulate the equations that describe the pair 
wind from a hot bare strange star, and the bound- 
ary conditions. In § 3 we describe the computa- 
tional method used to solve these equations. In 
§ 4 we present the results of our numerical simula- 
tions. Finally, in § 5, we discuss some astrophysi- 
cal applications. 

2. Formulation of the problem 

We consider, then, an pair wind that flows 
away from a hot, bare, unmagnetized, strange star 
with a radius of i? = 10^ cm. We assume that the 
temperature is the same over the surface of the 
star. The pair flux Fg from the unit surface of 
strange quark matter (SQM) depends on the sur- 
face temperature Ts only (Usov 1998, 2001a); and 
the pair wind is spherical. At a distance r from 
the stellar center and at time t, the state of the 
plasma may be described by the distribution func- 
tions fe±{p,i",t) and fj{p,r,t) for positrons (-I-), 
electrons (— ), and photons, respectively, where p 
is the momentum of particles. There is no emis- 
sion of nuclei from the stellar surface, and there- 
fore, the distribution functions of positrons and 
electrons are identical, /e+ = /e- = /e/2. 
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In the first run of our simulations tlie outflow 
of pairs from the stellar surface starts at the 
moment t = 0, and its rate is constant with time 
(E = lO^** ergs s~^) till the moment when a steady 
state is established in the examined space domain 
(see below). Then, the rate of the pair outflow in- 
creases instantly by a factor of 10, and the second 
run starts and continues till a new steady state is 
established, and so on. So, we find the structure 
of stationary pair winds and their emerging ra- 
diation for E = lO^^^, 10^5, ... , 1042 ergs s'^. 

2.1. Equations 

We use the relativistic Boltzmann equations 
for the pairs and photons, whereby the dis- 
tribution function for the particles of type i, 
/i(|p|, /i, r, i), (i = e for pairs and i = j for 
photons), satisfies (e.g., de Groot, van Leeuwen, 
& van Weert 1980; Mihalas 1984; Mezzacappa & 
Bruenn 1993) 

(1) 

where fj, = cos ^, ^ is the angle between the radius- 
vector r from the stellar center and the particle 
momentum p, p= |p|, and 

Pe=Ve/c=[l-{meC'/eef]^/^ P-y = 1 , (2) 

while 

Ce = c[p^ + (mec)2]i/2 and = pc (3) 

are the energy of electrons and photons, respec- 
tively. Also jj' is the emission coefficient for the 
production of a particle of type i via the physical 
process labeled by q, and xt is the correspond- 
ing absorption coefficient. The summation runs 
over all considered physical processes that involve 
a particle of type i. Gravity is neglected (but its 
effects are briefly discussed in §5). This is a good 
approximation for E > iEdd, while for E < ^Edd 
the gravity effects may be important, where 

LEdd ^ 7 X 10^\M/Mq) ergs s'^ , (4) 

is the Eddington luminosity for the pair plasma, 
above which radiation-pressure forces dominate 
over gravity. 



The particle densities are given by 

ni{r,t) ^ J Mp,r,t)dp. (5) 

For convenience of numerical simulations we 
use, instead of /», the quantities 

Ei{ei, iJ.,r,t) = ^3 , (6) 

as are used in the "conservative" numerical 
method, which can provide exact conservation of 
energy on a finite computational grid (see below). 
Since 

Cifidrdp = ' ' ' drdcidfi = Eidrdcidfi , (7) 

we see that Ei is the energy density in the 
{r,/i, ej,} phase space, in which the volume el- 
ement is drdcidii — Airr^drdeidiJ. 

Prom equations (1) and (6) the Boltzmann 
equations can be written in terms of Ei as 

+ ^^ir'P.E,) + - # [(1 - (^'^E.] 
c at r-^ or r on -' 

= ^[7??-X?S,], (8) 

Q 

where 

This form of the Boltzmann equations is said to 
be "conservative" , in the case without gravity, be- 
cause the "transport" terms can be written as a 
derivative of a flux divided by a volume (Harleston 
& Holcomb 1991). This allows one to develop 
a conservative finite difference code for numeri- 
cal simulations and to carry out calculations with 
large time steps even for very large optical depths. 
In our simulations we use the Boltzmann equations 
(8) for both pairs and photons. 

2.2. Boundary conditions 

The computational boundaries are 

r-int < r- < Text, < i < , (10) 

where Text = 1-65 x 10^ cm is chosen as the radius 
of the external boundary, and t^t is the time when 
the wind approaches stationary closely enough. 
The internal radius is that typical of a strange 
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star (rint = R)- At this radius, the injection rate 
of pairs is taken as (Usov 1998, 2001a) 



where 



3 X 10*^ exp (-0.593C) 



(11) 



ln(l + 2C' 

(1 + 0.0740^ ' 2(13.9 + 



+ 



cm ^s ^, 



(12) 



and C = 20(Ts/109K)-i. The energy spectrum 
of injected pairs is thermal with temperature Ts, 
and their angular distribution is isotropic for < 
M < 1- 

The energy injection rate is equal to the ther- 
mal luminosity of the strange star in pairs, 



E = Lss = N^mec^ + (3/2)fcBTs] . 



(13) 



where fee is the Boltzmann constant (Usov 2GGla). 
For the range of the energy injection rates we 
consider, E = 10^^ — 10^^ ergs s~^, the surface 
temperature is in a rather narrow range, Ts ~ 
4 X IQS - 109 K (see Fig. 1). 

The thermal emission of photons from the sur- 
face of a bare strange star is strongly suppressed 
if the surface temperature is not very high, Tg 4C 
10^1 K (Alcock et al. 1986a). The reason is that 
the plasma frequency of quarks in SQM is very 
large, fkVp ~ 20 — 25 MeV, and only hard pho- 
tons with energies > hiOp can propagate in 
SQM. The luminosity in such hard photons, which 
are in thermodynamic equilibrium with quarks, 
decreases very fast for Ts ^ hwp/kQ (Chmaj, 
Haensel, & Slomihski 1991; Usov 2001a), and in 
our case when Tg is < 10^ K, it is negligible (see 
Fig. 1). However, low-energy (e-y < fuvp) photons 
may leave SQM if they are produced by a nonequi- 
librium process in the surface layer of thickness of 
~ c/ujp ~ 10-^2 cm (Chmaj et al. 1991). Re- 
cently, the emissivity of SQM in noncquilibrium 
quark-quark bremsstrahlung radiation has been 
estimated by Cheng & Harko (2003), taking into 
account the Landau-Pomeranchuk-Migdal effect, 
and it is shown that the SQM emission in noncqui- 
librium photons is suppressed at least by a fac- 
tor of 10^ in comparison with blackbody emission 
(see Fig. 1). It was argued that SQM is a color 
superconductor with a very high critical temper- 
ature Tc ~ 10^^ - 10^2 K (Alford, Rajagopal, & 



Wilczek 1999; Alford, Bowers, & Rajagopal 2001). 
In this case, the noncquilibrium emission is sup- 
pressed further. In our simulations, we neglect the 
surface emission of photons altogether. 

The stellar surface is assumed to be a perfect 
mirror for both pairs and photons. At the ex- 
ternal boundary (r = Text), the pairs and photons 
escape freely from the studied region, i.e., the in- 
ward (/X < 0) flux of both pairs and photons 
vanishes. 

2.3. Physical processes in the pair plasma 

Although the pair plasma ejected from the 
strange-star surface contains no radiation, as the 
plasma moves outwards photons are produced by 
pair annihilation (e+e~ 77). Other two-body 
processes that occur in the outflowing plasma are 
M0ller (e+e"*" — » e+e"*" and e~e~ — * e~e~) and 
Bhaba (e+e~ — > e+e~) scattering, Compton scat- 
tering (76 7e), and photon-photon pair produc- 
tion (77 e+e~). Two-body processes do not 
change the total number of particles in a system, 
and thus c;annot, in themselves, load to thermal 
equilibrium. For this reason, we include radiative 
processes (bremsstrahlung, double Compton scat- 
tering, and three photon annihilation with their 
inverse processes) that change the particle num- 
ber even though their cross-sections are at least 
^ ^ 10^ times smaller than those of the 

two-body processes (a = e^/Hc = 1/137 is the 
fine structure constant). Radiative pair produc- 
tion, 77 ^ e'^e^7, which is a radiative variant 
of photon-photon pair production, 77 e+e", 
is not essential and is ignored in our simulations. 
The processes we include are listed in Table 1. 

3. The computational method 

In the numerical scheme we define a grid in 
the {r, /i, e} phase-space as follows. The r domain 
{R< r < Text) is divided into jmax spherical shells 
whose boundaries are designated with half integer 
indices. The j shell (1 < j < jmax) is between 
rj-1/2 and rj+1/2, with Ar^ = rj+1/2 - rj-1/2 



R and r. 



.-1-1/2 



^-ext)- 



The /i-grid is made of fcmax intervals A/Ufc = 

- Mfc-1/2: 1 < < fcmax- 

The energy grids for photons and electrons are 
different. They are both made of uimax energy in- 
tervals Aci^oj = ei,u+l/2 - %u-l/2- 1 < W < Wmax, 
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but the lowest energy for photons is 0, while that 
for pairs is nieC^- 

The quantities we compute are the energy den- 
sities averaged over phase-space cells 

Ei,^,k,j{t) = ^ [ Eidediir^dr. (14) 

where AX = Aea;A/ifcA(r|)/3 and A{rf) = 

'j+1/2 'j-l/2- 

Replacing the space and angle derivatives in the 
Boltzmann equations (6) by finite differences, we 
have the following set of ordinary differential equa- 
tions (ODE) for Ei^a>,k,j specified on the computa- 
tional grid: 

1 dEi^cj^kj A[r'^fj,kEi^^^k]j 



c dt 



+ 



A(rf)/3 
1\ ^ A[{1-^^)E^^ 



(15) 



where = 1, 

Pe,u = s/l^^^^^^^^^/Z^ , (16) 

ei,u; = (ei,ij-i/2 + ei,w+i/2)/2 , (17) 

Mfc = (Pk+1/2 + Mfe-i/2)/2 , (18) 

1\ (^1+1/2 - r-|_l/2)/2 ^^g^ 



.r/j (rf+i/2-^|-i/2)/3' 

-r|_i/2[/^fc-2^i,a;,fe]r=rj_i/2' (^0) 

A [(1 - fJ''^)Ei,u,,j]k = (1 - fJ'k+l/2)[Ei,aj,j]ix=tJ.k+i/2 

~(1 ~ Mfe-l/2)[-E'i,w,j]A»=Mfe-i/2 > (21) 

Ei,u,k{r) = -r-^ — / Ei{€,iJ.,r)dedti, 



Ei,u;,j{l^) = . \ 3 / Ei{e,iJ,r)der'^dr 



(22) 



(23) 



[Mfe-E'i,w,fe]r=rj+i/2 — (1 ^ Xi,w,fcJ+l/2) 



Ei (jj k,j ~i~ 



+Xi,w,fc,j+l/2Mfc 



J ■ 2 



1 



, (24) 



,(25) 



+ 



[-^j,w,i]Ai=/jfc+i/2 ~ Ei i^ kj 
AiJ,k{Ei^^^k,j — £'i,aj,fc-i,j) 



A/i/s_i -I- A/x/s 



(26) 



'i,uj,k 



:> AX 



rfidediir^dr, (27) 



XiE^dedjir^dr, 



(28) 

For the physical processes included in our simu- 
lations the expressions '']'i^ui,k,3 ^'^'^ ixE)i^^^k,j 
given in the Appendix. 

The dimensionless coefficient x is introduced 

to describe correctly both the optically thin and 
optically thick computational cells by means of a 
compromise between the high order method and 
the monotonic transport scheme (see Richtmcycr 
& Morton 1967; Mezzacappa & Bruenn 1993; Ak- 
senov 1998). 

Energy conservation on the grid may be written 

as 



j2 ^dE.^^.kd A... 



(29) 



where 



Ej + l/2 — '^1+1/2 X/ ^^ul^^J'kPi,ul[^^kEi^u,k]r=rj+l/2 
i,uj,k 

is the energy flux through the sphere of radius 
rj+i/2. In our scheme this is exactly satisfied. Par- 
ticle number is also explicitly conserved when only 
two-body processes are taken into account; par- 
ticle production occurs only via three-body pro- 
cesses. 

There are several characteristic times in the sys- 
tem. Some are related to particle reaction times, 
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some to the time to reach steady state. These 
times may greatly differ from each other, espe- 
cially at high luminosities (> 10^^ ergs s~^), when 
the pair wind is optically thick. The set of equa- 
tions (15), which describes the optically thick wind 
is said to be stiff. (A set of differential equations 
is called stiff if at least some eigenvalues of Ja- 
cobi matrix differ significantly, and the real parts 
of eigenvalues are negative.) In contrast to Mezza- 
cappa & Bruenn (1993) we use the Gear's method 
(Hall & Watt 1976) to solve the Boltzmann equa- 
tions (15). This high order implicit method was 
developed especially to find a solution of stiff sets 
of ODE. To solve a system of linear algebraic 
equations at any time step of the Gear's method 
we use the cyclic rcdiiction method (Mczzacappa 
& Bruenn 1993). The number of operations per 
time step is oc (ct;maxfcmax)^imax, which increases 
rapidly with the increase of ciJniax and fcmax • There- 
fore, the numbers of energy and angle intervals are 
rather limited in our simulations. 

Below, we use the e, /x, r-grid with Wmax = 
fcmax = 8 and jmax = 100. The discrete energies 
of the e-grid are 

2 meC^ f 1 - COSTTCCo;, X^i < I, 

'''^ ^'^ ' 2 \ 3 -h cosTraJo;, x^, > 1, 

(30) 

and ei,a;„3,^+i/2 = oo, where = (w - l)/(u;max - 
3). This gives a denser grid at low energies and 
near the threshold of pair production, e = rricC^. 
For photons (m^ = 0) the discrete energies (in 
keV) of the e-grid are 0, 48.8, 176.5, 334.4, 462.1, 
510.9, 559.7, 687.4, and oo, i.e., three energy in- 
tervals are above the pair production threshold. 
The /i-grid is uniform, A/Xfc = 2/fcinax = 1/4. 
The shell thicknesses are geometrically spaced, 
Arj = 1.3Arj_i, and the thickness of the initial 
shell is An = 2 x 10~^ cm. 

We ran two test problems. In the first we as- 
sumed that the external boundary is a perfect mir- 
ror for both pairs and photons. Pair injection 
was stopped at some moment, and the subsequent 
evolution of the system was followed. We found 
that it evolved eventually to the state when the 
distributions of photons and pairs are completely 
isotropic and uniform. In the second we verified 
that the processes we included lead to thermal 
equilibrium in a spatially uniform system (this also 
checks the adequacy of our energy grid). Follow- 



ing Pilla & Shaham (1997), we considered the time 
evolution of photons and pairs that start far from 
thermal equilibrium but have isotropic and uni- 
form distributions. At f = the total energy den- 
sity of photons and pairs was taken to be that 
in thermal equilibrium at T = 10'' K {U-qT ~ 100 
keV) . The initial photon and the pair distributions 
were flat and nonzero within the the energy inter- 
val 10~^ < (cj - mi(?)lme.(? < 2. The initial den- 
sities of photons and pairs were equal. As seen in 
Figures 2 and 3, at time t = 2.7x 10~* s, when the 
system is almost stationary, photons and pairs are 
near thermal equilibrium. For photons and elec- 
trons, the differences between the calculated and 
equilibrium spectral energy densities are a few tens 
perccnts (attributed to the coarseness of the grid) 
which is then the accuracy of calculations of en- 
ergy spectra in our simulations. 

To check the effects of grid coarseness we also 
performed test computations with fcmax = 4 and 
6, and separately with Wmax = 6 and 7 for the 
discrete energies of the e-grid given by equation 
(30) where x^, = {oj - l)/(wmax - 2), i.e., only 
two high-energy intervals are above the pair pro- 
duction threshold. We did not observe essential 
changes in the results. This implies that the main 
results of our simulations are not sensitive to the 
e, yu-grid in spite of its being rather coarse. 

4. Numerical results 

We give here the results for the properties of 
pair winds. The energy injection rate, E, is the 

only parameter we vary in our simulations. As 
explained in § 2, we start from an empty wind 
injecting pairs at a rate 10^^ ergs s~^. After a 
steady state is reached we start a new run with 
this steady state as initial condition, increase the 
energy injection rate by a factor of 10, wait for 
steady state, and so on. Figure 4 shows the to- 
tal emerging luminosity in photons and in pairs at 
the external boundary as a function of time t. We 
see that in each run this luminosity increases even- 
tually to its maximum value L, which equals the 
energy injection rate, L — Le + = E. The rise 
time is ^ 10""^ — 10~^ s, which is the characteristic 
time on which the pair wind becomes stationary in 
the examined space domain. We next present the 
results for the structure of the stationary winds 
and their emerging emission. Figure 5 shows the 
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mean optical depth for photons, from r to r^xt, 
defined as 



where 



i: 



dr 



J de d^ x-y (e, 



/ dt djd E^{e, fi, i 



X7 



E 



X-y ' 



(31) 



(32) 



The contribution from Text to infinity is negligi- 
ble for r < 10® cm, so r^(r) is practically the 
mean optical depth from r to infinity for these 
values of r. The pair wind is optically thick 
[t^{R) > 1] for E > 10^7 ergs s^^. The radius 
of the wind photosphere rph, determined by con- 
dition T(rph) = 1, varies from ^ R for E = 10^'' 
ergs s~^ to ~ lOR ~ 10'' cm for E = 10^^ ergs s~^. 
The wind photosphere is always deep inside our 
chosen external boundary (rph < O.lrcxt), justi- 
fying our neglect of the inward (/x < 0) fluxes at 

r = Text- 

The emerging luminosities in pairs {Le) and 
photons (Z/-y) are shown as fractions of the total 
luminosity L = + in Figure 6. For L < 
icq — 2 X 10^^ ergs s~^, the injected pairs remain 
intact, and they dominate in the emerging emis- 
sion (Z/e > Lj). For L > Leq, the emerging emis- 
sion consists mostly of photons {Lj > Lg). This 
simply reflects the fact that in this case the pair 
annihilation time tann ~ {neO'Tc)~^ is less than 
the escape time tesc ~ R/c, so most injected pairs 
annihilate before they escape (ctt is the Thomson 
cross section). The condition fann — iesc implies 
(e.g., Beloborodov 1999 and references therein) 

L ~ ^![!!^ 2.3 X 10^5 ergs s'^ , (33) 

(Tt 

which practically coincides with Lgq- For very 
high luminosities (L ^ L^q), reconversion of pho- 
tons into pairs is inefficient, as the mean photon 
energy at the wind photosphere is rather below 
the pair-creation threshold (see below), and pho- 
tons strongly dominate in the emerging emission. 

Figure 7 shows the rates of outflow of pairs 
(A^'e) and photons (N^) through the surface at ra- 
dius r as functions of r. At r = R, the pair outflow 
rate is put equal to the rate of pair injection 



while the photon outflow rate is zero (iV^" = 0). 
(Because the surface temperature depends weakly 
on E, N™ is nearly proportional to E.) There 
is the upper limit on the rate of emerging pairs 
jy-max ^ ^qAS g-i_ jf > TV^max^ ^^le rate of pafr 

outflow Ne decreases sharply at the distance 



^ann — 



1 



~ 10" 



E 



cm 



lQi2 gj-gg g-1 

. . ' (35) 

from the stellar surface, where ^ = N™ / AttR^v°^^ 
is the density of pairs at the surface, v°^^ ~ 10^ 
cm s~^ is the velocity of the pair plasma outflow 
near the surface (see below), (Tann — '^t{c--/ {ve)) 
is the pair annihilation cross section, and {ve) — 
(3fcBTs/me)^/^ ~ 10^° cm s~^ is the mean velocity 
of injected pairs. As a result, the rate of emerging 
pairs is limited to N^'^^ within a factor of 2. For 
E > 10^* ergs s~^, radiative three-body processes 
arc important, and the total rate of the parti- 
cle outflow N = Ne + increases with radius 
(see Fig. 7). For E = 10^^ ergs s~^, when pho- 
tons strongly prevail in the emerging emission, the 
rate of emerging photons increases by a factor of 
15 in comparison with N™. The rates of energy 
outflow in pairs (Eg) and photons (E^) vary 
with radius more or less similarly to the particle 
outflow rates, except that the total energy rate 
(Eg + E^) doesn't depend on radius and equals to 
the energy injection rate E (see Fig. 8). The up- 
per limit on the emerging luminosity in pairs 
is L™'^'^ ~ 2 X 10'^'' ergs s~^ (which includes the 
rest energies). The existence of the upper limits 
^max g^j^j L™^^ are connected with production of 
photons via radiative three-body processes. This 
can be seen in Figure 9, which shows the rates of 
particle outflow when only two-body processes are 
included. We see that for E 10'^^ ergs s""'^ and 
r — R^ /ann the outflow rate of pairs is only 
a few times less than that of photons. 

The bulk velocity of the pair plasma outflow 



(36) 



7V^° ~ 10^«(E/10^2grgg g-i) g- 



(34) 



is shown in Figure 10, where rig is the pair number 
density. We see that at low luminosities (E < 10^'' 
ergs s~^) increases with distance from the sur- 
face, reaching, for E = 10'^'' ergs s~^, as high a 
value as 0.7c. For such luminosities the wind is op- 
tically thin, r-y < 1 (see Fig. 5), and pairs and pho- 
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tons flow away more or less independently. The in- 
crease of v°^^ occurs because outflowing pairs are 
heated by annihilation photons via Compton scat- 
tering. At high luminosities (E > lO''^ ergs s"'^) 
the velocity of the pair plasma outflow decreases 
at r — i? > lann and is as small as ~ 2 x 10^ cm 
for E ~ 10*1 - 10^^ ergs s'^ (see Fig. 10). To ex- 
plain the decrease of we introduce the bulk 
velocity of the photon gas outflow 



(37) 



where is the photon number density. The ra- 
tio v°'^^/c is a measure of photon anisotropy and 
varies from zero at the stellar surface to 1 far from 
the wind photosphere (see Fig. 11). In our case 
when the injected plasma consists of pairs with 
mean velocity (ve) ~ c, the free path length of 
Compton scattering, which is the main mecha- 
nism of opacity for photons, is of the order of the 
free path length of pair annihilation Zann- There- 
fore, at the distance of ^ l^nn from the stellar 
surface the outflowing pair plasma is decelerated 
by nearly isotropic {v°^^/c <C 1) photons. Then, 
the pair plasma, which is optically thick, flows 
away with nearly the same velocity as the pho- 
ton gas (?;°"* ~ u°"*) up to the wind photosphere 
(see Figs. 10 and 11). It is worth noting that for 
E > 10^^ ergs s"^ and r — i? » Zann, when photons 
prevail, the wind dynamics is mostly determined 
by photons, pairs being only responsible for pho- 
ton opacity. 

Figure 12 shows the pair number density as a 
function of the distance from the stellar surface. 
As seen in Figures 7, 10 and 12, for iVf > iV™^'' 
the rate of pair number outflow = 47rr^neW°"* 
decreases very fast with radius at r — R > /ann 
because both Ue and decrease fast at the same 
distances. 

Figures 13 and 14 show, respectively, the mean 
energies of photons and pairs, which are deter- 
mined by 



J eifi{€i,iJ.,r) dcidfi 



(38) 



as functions of the distance from the surface. We 
see that the mean photon energy (e-^) decreases 
with the increase of E at all radii. In the E range 
from ~ 10^* to ^ 10^'^ ergs s^^, where most of 
the photons in the system are produced via pair 



annihilation, the decrease of (e-y) is rather weak 
and occurs because of energy transfer from anni- 
hilation photons to pairs via Compton scatter- 
ing. As a result, the emerging pairs are heated 
up to a mean kinetic energy (ce) — 500 keV at 
E ~ 10^'^ ergs s'^ (see Fig. 14). For E > 10^^ 
ergs s~^, (e-y) decreases mainly because of produc- 
tion of rather low-energy photons in the radiative 
three-body processes (see Fig. 15). For E = 10*^ 
ergs s~i we have (e-^) ^ 40 keV for the emerg- 
ing photons. This value of (e-y) is near the mean 
energy of blackbody photons for the same energy 
density as that of the photons at the wind photo- 



sphere, which is (e. 



bb\ 
7 ' 
,bb\ 



30 keV. (The difference 



between (e-y) and {e^ ) is less than the energy res- 
olution of the e-grid at low energies, which is 20 
keV.) 

Our Figure 15 shows that for high luminosities 
{L > 2 X 10'^^ ergs s""^) the mean energies of pho- 
tons and pairs decrease with luminosity. This is 
consistent with previous studies showing a similar 
behavior for large enough values of the compact- 
ness parameter 



2tt— (39) 



(Svensson 1984 and references therein). 

Figures 16 and 17 present the energy spectra of 
the emergent photons and pairs. At low luminosi- 
ties, E ~ 10^'' — 10'^^ ergs s"^, photons that form 
by pair annihilation escape more or less freely, and 
the photon spectra resemble a very wide annihila- 
tion line. For E > 10^* ergs s"^, changes in the 
particle number due to radiative three-body pro- 
cesses arc essential, and their role in thcrmaliza- 
tion of the outflowing plasma increases with the in- 
crease of E. We see in Figure 14 that for E = 10*^ 
ergs s""'^ the photon spectrum is near blackbody, 
except for the presence of a high-energy tail at 

> 100 keV. Such a hard spectrum of photons 
together with the anti-correlation between spec- 
tral hardness and photon luminosity (see Fig. 15) 
could be a good observational signature of a hot, 
bare, strange star. 

The high-energy tail of the emergent photons 
covers six e-grid intervals and is real. Moreover, 
if the energy injection rate E is less than 10*^ 
ergs s~i, such a high-energy tail has to be present, 
i.e., the energy spectrum of the emergent photons 
cannot be completely a Planckian at high energies. 



8 



Indeed, if at the wind photosphere, photons and 
pairs are in thermal equiUbrium, their tempera- 
ture is 

For E = 10"*^ ergs s"""^, when rph is ~ 10^ cm (see 
Fig. 5), we have T^^ ~ 0.61 x 10^ K. At this tem- 
perature the density of equiUbrium pairs is (e.g., 
Paczynski 1986) 

(41) 

or numerically, rig.ph ^ 10^^ cm while for 
E > 10^^ ergs s^^ the pair density at the pho- 
tosphere cannot be essentially smaller than ^ 
l/(rphCTa„n) ~ 10^^ - IQis cm-3 (see Fig. 12), 
at which the optical depth for pair annihilation 
is ~ 1. Even if we take the stellar radius as 
the radius of the wind photosphere, rph = R, 
for E = 10^^ ergs s~^ from equations (40) and 
(41) we have Tp^J ~ 1.93 x 10^ K and ne,ph ^ 
10"'^^ ^ l/(rphO'ann)- Hcncc, in this case the en- 
ergy spectrum of the emergent photons cannot be 
a Planckian at high energies, and it has to range 
to ^ iTieCp. This spectrum might be completely a 
Planckian starting only from E ~ 10'''^ ergs s^^, 
at which for r ~ i?, we have ~ 3 x 10^ K and 

> l/(rphf7ann)- 

5. Discussion 

We have identified certain characteristics of the 
expected radiation from hot, bare, strange stars 
that, we hope, will help identify such stars if they 
exist. The spectrum, we find, is rather hard for 
the studied luminosity range. This makes such 
stars amenable to detection and study by sensi- 
tive, high energy instruments, such as INTEGRAL 
(e.g., Schoenfelder 2001), which is more sensitive 
in this range than previous detectors. 

As the emission from bare, strange stars 
is characterized by super-Eddington luminosi- 
ties (see Fig. 1) and, as we now find, also by 
hard X-ray spectra, soft 7-ray repeaters (SGRs), 
which are the sources of short bursts of hard X- 
rays with super-Eddington luminosities (up to 
~ 10^^ — 10**^ ergs s~^), are potential candidates 
for strange stars (e.g., Alcock, Farhi, & Olinto 



1986b; Cheng & Dai 1998; Usov 2001b). The 
bursting activity of SGRs may be explained by fast 
heating of the stellar surface up to the temperature 
of ^ (1 — 2) X 10^ K and its subsequent thermal 
emission (Usov 2001b, c). The heating mechanism 
may be either impacts of comets onto bare strange 
stars (Zhang, Xu, & Qiao 2000; Usov 2001b) 
or fast decay of superstrong (~ 10^^ — 10^^ G) 
magnetic fields (Usov 1984; Thompson & Duncan 
1995; Heyl & Kulkarni 1998). For typical lumi- 
nosities of SGRs (L ~ 10-*! - lO'^^ ergs s'^), the 
mean photons energy we find is ~ 40 keV (see 
Fig. 15), which is consistent with observations of 
SGRs (Hurley 2000). The rise time (~ lO'^-lO'^ 
s) of the emerging luminosity (see Fig. 4) implies 
that the variability of bursts from SGRs may be 
explained in the strange star model. 

Another important idiosyncrasy that we find is 
a strong anti-correlation between spectral hard- 
ness and luminosity. While at very high luminosi- 
ties {L > 10^^ — 10^^ ergs s~^) the spectral tem- 
perature increases with luminosity as in blackbody 
radiation, in the range of luminosities we studied, 
where thermal equilibrium is not achieved, the ex- 
pected correlation is opposite (see Fig. 15). Such 
anti-correlations were, indeed, observed for SGR 
1806-20 and SGR 1900-1-14 where the burst statis- 
tic is high enough (e.g., Feroci et al. 2001; Gogus 
et al. 2001; Ibrahim et al. 2001). This is encour- 
aging, but a direct comparison of this data with 
results such as ours will require a more detailed 
analysis. In particular, the effects of strong mag- 
netic fields have to be included. The processes 
we consider can be significantly modified by such 
fields, especially above 10^"^ G (e.g., Daugherty & 
Harding 1989). Qualitatively new processes such 
as photon splitting (e.g., Adlcr 1971; Baring & 
Harding 2001; Usov 2002) may become important 
in a strong magnetic field. The observed periodic 
variations in the SGRs may be due to a rotation 
of a star with non-uniform surface temperature, 
while our result apply for isotropic emission. We 
hope to deal with these effects elsewhere. 

We noted that there is an upper limit to the rate 
of emerging pairs N^'^^ ~ 10**^ s~^. Positrons can 
then annihilate in an ambient medium. Hence, 
the luminosity in annihilation emission from the 
region far from a hot bare strange star may be as 
high as 2mc^N™^ ~ 10^^ ergs s~-^. 

In our simulations, gravity is neglected. Grav- 
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ity decelerates the pair wind and increases the pair 
density near the surface. Also, photons are red 
shifted when they escape from the star's vicinity. 
To estimate the effects of gravity we ran simula- 
tions for a star mass of IAMq, where gravity was 
included in the Newtonian approximation. We 
found that at low luminosities (E < I/Edd — lO'"*''' 
ergs s~^) the pair density rig does indeed increase 
significantly near the surface, and for E = 10^^ 
ergs the value of rif, is about ten times higher 
than without gravity. At such low luminosities 
the velocity of the pair plasma outflow deceases 
by about a factor of 3 because of gravity effects 
and is equal to 10^ cm s~^ within a factor of 2 
(cf. Fig. 10). The probability of pair annihila- 
tion increases because of the pair density increase, 
and photons dominate the emerging emission for 
i = E > 3 X 10^4 ergs s'^ (cf. Fig. 6). At high 
luminosities (E > I/Edd) gravity doesn't affect 
the pair wind structure significantly, especially at 
E » I/Edd- For E < 10^'^ ergs when the pair 
wind is optically thin (see Fig. 5) the mean energy 
of emergent photons (cj) deceases by ~ 20% be- 
cause of the red shift, while for E > 10^^ ergs 
the decrease of (e-^) is less than a few pcrcents, 
compared with the accuracy of our simulations, 
which is not higher than ~ 10 — 20% because the 
e, /i-grid is rather coarse. Note also that gravity 
corrections in the Newtonian approximation are 
valid for the relativistic Boltzmann equations with 
an accuracy of a factor of two, and, while cap- 
turing some of the effects of gravity, they are not 
completely self-consistent. We have thus preferred 
to present our results for the self consistent case 
without gravity, deferring the full inclusion of rel- 
ativistic gravity for a later treatment. 

We are grateful to the anonymous referee for a 
careful reading of the manuscript and for many 
helpful comments. This work was supported 
by the Israel Science Foundation of the Israel 
Academy of Sciences and Humanities. 
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APPENDIX 



A. Emission and absorption coefficients for two- and three-body processes 

We use two sets of independent variables p. r. t. and e,fj,,r,t. We take a uniform distribution of particle 
density inside the grid volume Aey,a;A/XfcA(r|)/3: 

Ei{e,iJ.,r) Ei^u,,k,j . ^ ^ ^ ^ ^ ^ 

ni{e,iJ.,r) = = for ei,(^_i/2 < e < ei,a,+i/2, Mfe-1/2 < M < /"fe+1/2, rj_i/2 < r < rj+1/2, 

(Al) 

where €i^a> = £1,^-1/2 + ^ei,a;/2. (We suppress the r-dependence of functions.) 

A.l. Compton scattering 

The time evolution of the distribution functions of photons and pair particles due to Compton scattering, 
76 j'e', may be described by (Ochelkov et al 1979; Berestetskii et al. 1982) 



^ = /'dk'dpdp'«;k',p';k,p[/-y(k',t)/e(p',t)-/-y(k,t)/e(p,i)], (A2) 

/ 'ye — v^'q' J 

= /'rfkdk'rfp't«k',p';k,p[/-y(k',i)/e(p',i)-/-y(k,i)/e(p,t)], (A3) 



dt 
dt 

where 

c(5(e^ — Ee — e!l — , , \ivifi\ 



is the probability of the process, 

\Mf,\^ = 2'{nhrrlmy< 



2 9 



2 2 

h 


(r 


2 2 


2 2 
H : 


2 2 

u — m^c^ 




2 2 


u — rn^a 






2 2 


2 2 
M — m^c 






4 


s — rn^c^ 



(A5) 

is the square of the matrix element Mfi, r,. = e'^ /{meC^) is the classical electron radius, s = (p + 6)^ and 
u = (p — t'Y are invariants, I = (e-y/c)(l, e^) and p = (ee/c)(l, /3ee) are energy-momentum four vectors of 
photons and electrons, respectively, dp = dedoe^ (3^ / , and do = d^dcj). 

As an example, wc calculate the value {xE)Z^Z'k ° ■ start this calculation by considering the negative 
term in equation (A2) , which is responsible for the Compton absorption of photons; 



/a/,(k,^)^-'-- 

I dt 



7c— *7'c' 

Substituting equation (A4) into equation (A6) we can obtain 



J (ik'dpdp'u;k',p';k,p/7(k,t)/e(p,i) . (A6) 



7e->7'e' J-^cec^Cg 
/ _ £66^,(1 - /?be-b^) 



do'^dp f^{Kt)Up,t), (A7) 



dt 

^"ee(l-/?be-b;) + e^(l-byb;)' -e--e^.7--7' 



where 

4= XTm e; = ee + e^-e' (A8) 



11 



b, = p,/p, = p'Jp', K = (/3eebe + e,h^ - e;b;)/(/3'e'J, and |M/,|2 = \Mfi\y[d'{2nhn 
For photons, the absorption coefficient in the Boltzmann equations (1) is 

where drii = deidoiefl3ifi/c^ = deidoiEi/{2nei). 

From equations (A2) and (A9), wc can write the absorption coefficient for photon energy density Ej 
averaged over the e, /x-grid with zone numbers u and k as 

(x-g)rrr = , r = dn^dn^do' 7, . (AlO) 

Similar integrations can be performed for the other terms of equations (A2), (A3), and we have 

^^'-''^ - Ae„.AM. - Ae„.A^, j^^-e^^ ^^^^edo, ^g^^^^^, , (All) 



■^7 e 



fig € A/^^ 



deedfie'n2 



1 



Aee,wA;Ufc 



Aee,a;A/Xft 



dn-yduedo' — — , 
^ ^ 16eee^ 



(A12) 



(»e6An^ 



A A — = ~7 — / dn-yduedo'— — (A13) 

Aee,^Ank Aee,wA/Xfe yee€Aee,„ ' lee^-e; 

The emission and absorption coefficients in the Boltzmann equations (12) are given by 

<^,^.- = A^P <-^^'* ' = A(;f)75 (^^)-.^^''^^ ' ^^^4) 

where is 76 — > 7'e' for Compton scattering. 

To integrate equations (A10)-(A13) numerically over <j) {0 < (j) < 27r) we introduce a uniform grid ^;ipi/2 
with 1 < / < Imax and A(j)i = {(pi+1/2 — 4'i-i/2)/'^ = 27r//max- We assume that any function of ^ in equations 
(A10)-(A13) in the interval Adj is equal to its value ai (j) — (pj — {4>i-i/2 + 0j+i/2)/2- Since the problem 
is axi-symmetric it is necessary to integrate over (j) only once at the start of calculations. The number of 
intervals of the (/>-grid is taken as Zmax = 2fcinax = 16. 

A. 2. Two-photon pair annihilation and creation 

The rates of change of the distribution function due to e~e+ ^ 7172 are 

^ g " ) = dkjdp-dp+Wu,,]^2;P-,P+fe-{P-,t)fe+i.P+,t) , (A15) 



f^f22^11^\ = - f dkjdp-dp+wp_^p^.k^,kj^^{ki,t)f^{k2,t) , (A16) 



for i = 1, i = 2, and for j = 1, i = 2. 
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f dUip±,t) \ =- /■dp^dki(ik2t«ki,k.;p_,p+/e-(P-,i)/e+(P+,i), (A17) 

I TT 1 = Mp=Fc;kidk2w;p_,p+;ki,k2^^^ ir^ -, (A18) 

V / 7172— >e~e+ J 



where 



Wp_,p+;ki,k2 = ^ 77r^2 P-+P+-ki-k2 — — , (A19 

fe+ — fe- = /e/2 and /-^j 2 = f'y/2. Here, the matrix element is given by equation (A5) with the new 

invariants s = (p_ — 61)^ and u = (p_ — 62)^ (Berestetskii et al. 1982). 
The energies of photons created via annihilation of a pair are 

, m^c* + e_e+(l - /3_/3+b_-b+) , / » ^,^n 

^^(^^^= e-(ll/?-b-b") + e.(i%.b!.bO ' ^M=e-+e,-e,, (A20) 

while the energies of pair particles created by two photons are 

^_^^_^^BTV^HZ^^ e+(b_)=ei + e2-e_, (A21) 

where A = (ei + 62)^ - [(eibi + e2b2)-b_]2, B = {ei + 62)6162(1 - bi-ba), C = mlc^[{€ihi + e2b2)-b-]2 + 
^1^2(1 ~'3i''32)^- The sign in equation (A21) has to be chosen so that momentum is conserved in the reaction. 
Integration of equations (A15)-(A18) yields 

<J,"r"^"^= A \ if. dn,-dn^,do,§^+ f dn,-dn,,do,^^^] , (A22) 

7,u;,K Ae^,a,A/Xfc Wei6Ae^,a, 166-6+62 Je2£Ae^,„ 166_6+ / 



(X-E)e,a,,fe ^'^' = -r — / _ dn^-dn^+doi ' + / _ dn^-dn^+doi- 



A6e,a;AyUfc W^-eAc.,„ 166+62 J '+<^^'c.^ 166-62 



(A23) 



ixE)l?, r = 1 / dn^,dn-y^do-—^ — 11}-+ / dn^.dn^^do ' ' 



16616+ / 

(A24) 

/ e_/3_|M,,'^^ 

1661626+ 7t+eAee,u 



jyZVjr^^ ^ = — -7 — I / dn-y.dn-y^do- ' — ''— + I dny.dny.do- 



(A25) 

where drigT = c^Cqr^Oq^e^/Jqi/eT, and dn-y^ ^ = <^ei,2<^oi^2ei,2/7i,2- 

A. 3. M0ller scattering of electrons and positrons 

The time evolution of the distribution functions of electrons (or positrons) due to M0ller scattering, 
e^e^ — > e^e^, are described by 



dt 



^ ^= j '^Pj^^PWP2Wpi,p^;pi,P2[/ei(pi,^)/e2(P2>^) - /ei(Pl,i)/e2(P2,t)] , (A26) 
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with i = 1, j = 2, and with j = 1, i = 2, and where 

cS{ei +€2- e[ - e'2) 



Pi,P2;pi.P2 



(A27) 



\Mf,\' = 2%nhrrlmyi^ 



+ 4m^c^{g-mic^) 



+ 



+ Aml&'{u-mlc^) 



/ei(Pl,i) = /e(Pl,i)/2, /e,(p2,i) = /e(P2,t)/2, s = (pi + = 2{ml^ + Prp^) , 9 = (pi - Pi)' 

2(m^c^ - prpi), and u = (pi - P2)' = 2(m^c^ - prpj) (Berestetskii et al. 1982). 

The energies of final-state particles are 



,/ R2 _ An 

e[{h[) =BT , e'2{h[) = e, + - e[ 



(A29) 



where i = (ei + €2)' - (ei/3ibi-b; + e2/?2b2-b'i)2, B = (e^ + e2)[m2c^ + £162(1 - /?i/32brb2)], and C = 
m2c^(ei/3ibrbi + e2/32b2-b;)2 + [m^c^ + £162(1 - /3i/32brb2)]2. 

Intergration of equations (A26), similar to the case of Compton scattering, yields 



6162^6163 

'le,oj,k 



dne, dne2 do^ ——^ — 

1D£i£2£2 



ixE) 



e,uj,k 



dneidne^do-^ — — — ; 

16£2£2 



,^^^^Jn^,dn^Jo,^^-— 
, , , , e'i/3!|Alf»K 

dUei dries "Ol 77 7 

2eAee,u; 16£l£2 



(A30) 

(A31) 



where (inei_2 = rfei,2rfoi,2ei,2A,2/ei,2- 

A. 4. Bhaba scattering of electrons on positrons 

The time evolution of the distribution functions of electrons and positrons due to Bhaba scattering, 
e~e+ — > e~'e"'"', is described by 



' dfe±{P±,t) 
dt 



] = / dp=Fdp'_dp+«;p'_,p^;p_,p^[/e-(p'_,i)/e+(pV,i) - /e-(P-,i)/e+(P+,*)], 



(A32) 

where 



PL.P^;p-.P. ^2 ^^(P- + P+ - P- - pV) i6£:£;£^_£V ' (^33) 

\Mfi \ is given by equation (A28), but the invariants arc s = (p_ — p'^)^, 5 = (p+ — p'^)"^ and u = (p_ +p+)^. 
The final energies e'_ , e'_^_ are functions of the outgoing particle directions in a way similar to that in Section 
A.3 (see Berestetskii et al. 1982). 

Integration of equations (A32) yields 

=+^»-'e+' _ 1 f f e'^(3'_\Mf,\' , f e'_P'_m,\^\ 



vt , fh ^ = — / dtie- dne+ do_ —— -. h / , dUe- dne+ do'_ ^ „ , , 

A£e,^A/Xfc W=-^^-- 16£_£+£'+ y^'+eA.,,^ 16£_£+ 

(A34) 
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' ' Aee,a,A/ifc We-eAee,„ 16e+e^ Je+eA.,,„ 16e_e:^ / 

(A35) 

where drieT = d€^do^e?^(3zpfeT ■ 

A. 5. Radiative three-body processes 

Compton scattering and scattering of pairs can equalize the temperatures of photons, electrons, and 
positrons. In the absence of three-body processes, the total particle number is conserved, and the equilibrium 
distribution function of the photons need not take the Planck form. The equilibrium of the reaction e+e" 
7i72 leads to equality of the chemical potentials of pairs and photons, = ip^. 

When only two-body processes are included (p^ ^ 0, in general. To get more realistic spectra we should 
include reactions that don't conserve particle number. The rates of such three-body reactions are at least 
~ 10^ times smaller than the rates of the two-particle reactions. The three-body processes we include 
in our study are listed in Table 1. We adopt the following emission coefficients for these reactions (cf. Haug, 
1985; Svensson, 1984). 



1. Bremsstrahlung, ee — > ee'y 



4^(11.2 + 10.46'2) + - 



r?; « T = ne+ne-^ In 



5 J e^p{l/e)K2il/e)' 

(A36) 

8cr2 



4ai + 10.4.) (l + ^)] 2(V2 + 2. + 2.^,-5iEfc|5^, ,A37) 

where x = e/(TOeC^), C = e"" '''^^^, 9 = kBT^/ (171^0^), and K2{l/0) is the modified Bcssel function of 
the second kind of order 2. The electron temperature, Tg, for the energy density Eg may be determined 
from 

y Eednde = J El'^dnde , (A38) 



fe^i^) = tA^ ^^P (-'-T^ ) ■ (A39) 



where E^'i = 2'Ke^ffi/(?, 

^^"'^^ " (2^/i)3 V kT 

is the equilibrium distribution function for electrons and positrons (degeneracy of particles is neglected), 
and is the chemical potential. 

2. Double Compton scattering, je — > 767 

3. Three photon annihilation, e~e+ 777 

= exp(-x/^) . (A41) 

We use the absorption coefficient for three-body processes written as 

X^P = ^7^7^^;^ (A42) 

where 77^^ is the sum of the emission coefficients of photons in the three particle processes, E^"^ = 27re^/^*i/c^, 
and 

2 

= {2nhf[exp{e/kT)-l] ^^"^^^ 
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is the equilibrium distribution function for photons. 

Prom equation (8), the law of energy conservation in the three-body processes is 

[j2{v?-X?Ei)dt^de = 0. (A44) 

For exact conservation of energy in these processes we introduce the following coefficients of emission and 
absorption for electrons: 



3p ^ /« - x'^'' E^)dedi^ 



/ Eededjji 



or 



Ee J Egdedn 



ril^ = 0, if y (?7^P - x^^E^)dedn > , (A45) 
Xe' = 0, ifJ{v^^-Xj^E^)dedf,<0. (A46) 
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Table 1; Physical Processes Inehided in Siirmlations 



Basic Two-Body 


Radiative 


IntcTaetion 


Variant 


Moller and Bhaba 




scatt(;ring' 


Bremsstrahlung' 


(:(■ 


(-( ^ CC'i 


Compton scattering 


Double Compton scattering 


je 76 


76 767 


Pair annihilation 


Three photon annihilation 


e+e~ 77 


e~'~e~ •<-> 777 


Photon-photon 




pair production 




77 — > 6+6~ 
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Fig. 1. — The luminosities of a hot, bare, strange star in e+e^ pairs (dotted line), in thermal equilibrium 
photons (dashed line), and the total (solid line) as functions of the surface temperature T^. The upper limit 
on the luminosity in non-equilibrium photons, Lneq ^ IO^^Lbb, is shown by the dot-dashed line, Lbb being 
the blackbody luminosity. 
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200 400 600 800 

[keV] 

Fig. 2. — Time evolution of the photon spectrum in the test problem of a homogeneous distribution starting 
away from equilibrium. The dashed line is the black body spectrum with the temperature of 10^ K. Times 
are marked on the curves. 
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Fig. 4. — The total emerging luminosity in photons and pairs at the external boundary (r = Text) as a 
function of time in runs with different values of E (equal to the long-time asymptotic value of L). 
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Fig. 5. — The mean optical depth for photons, from r to infinity, as a function of the distance from the 
stellar surface, for different values of E, as marked on the curves. 
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Fig. 6. — The fractional emerging luminosities in pairs (Le) and photons (L-y) as functions of the total 
luminosity, L = Lg + Lj. 
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Fig. 7. — The rates of particle number outflow in photons (sohd hues), in e+e~ pairs (dashed hnes), and the 
total (dotted lines) through the surface at radius r as a function of the distance from the stellar surface for 
different values of E, as marked on the curves. 
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Fig. 8. — The rate of energy outflow, as in Figure 7. 
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Fig. 9. — The same as Figure 7, with only two-particle processes taken into account. 
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Fig. 10. — The velocity of the pair plasma outflow as a function of the distance from the stellar surface for 
different values of E, as marked on the curves. 
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Fig. 11. — The velocity of the photon gas outflow as a function of the distance from the stellar surface for 
different values of E, as marked on the curves {v1^^/c is a measure of photon anisotropy). 
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Fig. 12. — The pair number density as a function of the distance from the stellar surface for different values 
of E, as marked on the curves. 
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g. 13. — The mean photon energy as a function of the distance from the stellar surface for different values 
E, as marked on the curves. 
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Fig. 14. — The mean kinetic energy of electrons and positrons as a function of the distance from the stellar 
surface for different values of E, as marked on the curves. 
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Fig. 15. — The mean energy of the emerging photons (thick solid line) and electrons (thin solid line) as a 
function of the total luminosity. For comparison, wc show as the dotted line the mean energy of blackbody 
photons for the same energy density as that of the photons at the photosphere. Also shown as the dashed 
line is the mean energies of the emerging photons in the case when only two particle processes are taken into 
account. 
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Fig. 16. — The energy spectrum of emerging photons for different values of E, as marked on the curves. The 
dashed hne is the spectrum of blackbody emission. 
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Fig. 17. — The energy spectrum of emerging pairs for different values of E, as marked on the curves. For 

comparison, the dashed Hnc is the energy spectrum of pairs that move with the velocity 77°"* ~ 4.6 x 10^ 
cm s~^ and have a Maxwellian spectrum with a temperature of 1.7 x 10^ K in the pair plasma frame. 
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